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Abstract 

We investigate an optomechanical system in which a cavity with a moving mirror is driven by two 
external fields. When the field frequencies match resonance conditions, we show that there exists 
a class of dark states of the moving mirror in the single-photon strong-coupling regime. These 
dark states, which cause the cavity to be decoupled from the external fields, is a manifestation of 
quantum coherence associated with the mirror's mechanical degrees of freedom. We discuss the 
properties of the dark states and indicate how they can be generated by optical pumping due to 
the decay of cavity field. 

PACS numbers: 42.50.-p, 42.50.Gy, 42.50.Wk 
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I. INTRODUCTION 



Quantum effects of a moving mirror interacting with electromagnetic fields in a cavity 
via radiation pressure has been a subject of considerable research interest recently l|-|5]. 
The motivation of this research is not only because the optomechanical coupling can pro- 
vide novel applications in quantum information processing, such as the storage of optical 
information as a .necnanical excitation Q and opton— transducers f 0r long-distance 
quantum communication but also because the system could be a platform to explore 
fundamental quantum phenomena at macroscopic scales. These include, for example, quan- 
tum entanglement 1^], Schrodinger cat states 14-17], and the modification of uncertainty 



relations due to quantum gravity 



In this paper we show how a harmonically bounded mirror can evolve into a dark state 
when the cavity is driven by two laser fields at certain resonance frequencies. Dark states are 
well known in a A— type three-level atom, in which a coherent superposition of two atomic 
ground levels can suppress light absorption completely. Our dark states reported in this 
paper share a similar mechanism, i.e., by quantum interference an optomechanical cavity 
cannot absorb photons from the external lasers. 

We note that interference effects of an optomechanical cavity driven by two light fields 



have been discussed in literature 



19|-|24j. By linearizing the system equations, an optome- 



chanically induced transparency (OIT) effect, which is an analogy of electromagnetically 



induced transparency, 
served in experiments 



ras been studied theoretically 19l-l21|. Recently OIT has been ob- 



22 



24j. Different from previous theoretical work, here we focus on 



the single-photon strong-coupling regime in which the radiation pressure from a single pho- 
ton can displace the mirror with a displacement comparable to the zero-point fluctuations. 
In such a regime, the usual linearized photon-phonon theory becomes inadequate because 
of the significant quantum fluctuations of the fields and the mirror. To study the physics 
in the strong coupling regime, one needs to solve the full quantum dynamics beyond the 
linear approximation, and some authors have reported interesting features, such as photon 



blockade effect 
cooling [28 



251 ] . multiple mechanical sidebands [26|, single-photon scattering [27| and 



It should be noted that in previous studies of OIT, the transparency refers to the non- 
absorption of a probe field, while the cavity contains a large number of photons due to the 
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presence of a control field [1914241] . Here we show that by exploiting the single-photon strong 
coupling, the driven cavity can have zero photon when the mirror is in the dark state. In 
this paper we discuss how such dark states of the mirror exist under a certain rotating wave 
approximation, and we indicate how they can be prepared by an optical pumping effect. 

II. THE MODEL 

We consider an optomechanical cavity formed by a harmonically bounded movable end- 
mirror and a fixed end- mirror [Fig. 1(a)], in which the cavity field and the movable end-mirror 
are coupled with each other via radiation pressure. The optomechanical cavity is driven by 
two lasers with frequencies U\ and o^- The Hamiltonian of the system is given by 

H = u c a?a + u M b ] b - ga^a (6 f + b) + [(f^e^ 1 * + fi 2 e~^ 2 *) a f + h.c] (1) 

where a (b) and u c (ojm) are respectively the annihilation operator and resonant frequency 
of the cavity field (mechanical) modes. The fii and Q2 are proportional to the amplitudes 
of the external fields. The parameter g is the single-photon coupling strength between the 
intracavity photon and the mirror caused by radiation pressure. In the frame rotating at 
the frequency u c , the transformed Hamiltonian is: 

H r = u u b% - ga ] a (b ] + b) + [(Oie"^ 1 ' + Q 2 e~ lA2t ) a ] + h.c] (2) 

where the detunings Ai = u\ — u c and A 2 = 0J2 — w c are defined. 

The first two terms of H r corresponds to the Hamiltonian H Q without driving, and it can 
be diagonalized by introducing a displaced oscillator basis as: 

H = u M b ] b - go) a (b ] + b) = ^ £ n, P IVVp) (Vv p |- (3) 

n,p 

Here the eigenvectors \if) n ,p) — \ n ) c ® D (ng/u>M) \p)m = \ n )c ® \P( n )) a/' with n (p) being the 
cavity photon (phonon) number. The \p{n)) M denotes the n-photon displaced Fock state 
of the mirror via the displacement operator D (ng/ujM) = exp[^ (tf — &)]. The energy 
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eigenvalues of Hq are e n)P = puiM ~ T^ff ' "which depend nonlinearly on photon number n, 
and linearly on phonon number p. 

By using the eigenbasis of H , H r becomes 

H r = £ ^P \^, P ) (V>n,p| + [ A S ( fiie ~' Alt + fi 2e" iA2i ) \lp n , p >) (Vn-l, P | + fc.C.] (4) 

n,p n,P,p' 
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FIG. 1: (Color online) (a) Schematic diagram of an optomechanical system consisting of a fixed 
end mirror and a movable end mirror with two driving fields, (b) The coupling scheme between 
energy levels of the optomechanical system (only zero- and one-photon states are shown). Here 
each laser is used to establish a set of resonant transitions, and \p) M = |i?(0)) M and \p) M = \p(l)) M 
for simplicity. By choosing g = gjy, there is no transition between |0) c |iV)M and \1) c \N)m- 

where we have expressed the annihilation operator a in the eigenbasis as 

a = Yl Y \^ n ^ ^"<pI a (VV,j/| = Y A pi' IV'n-l.p) (VVp'l ( 5 ) 

n,p n',p' n-iViV' 

with the coefficients A^ p/ = <fn (p\ £> f [(n - 1) g/u m ] D (ng/u M ) \p') = y/n (p\ D (g/u M ) \p')- 
Specifically, the coefficients are given by 

v^\/IM (-o p '~ p ^'- p (a ,p < p' 



An) 



where £ = g/u>M and L s r (x) are associated Laguerre polynomials. 



(6) 



III. EFFECTIVE RESONANT HAMILTON! AN WITH FINITE DIMENSIONS 



In this section, we show how evolution of the system state can be confined to a finite 
dimensional subspace involving only the zero and one cavity photon number and N displaced 
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phonon number states, assuming the initial state is the ground state. This is achieved by 
exploiting resonances and specific values of optomechanical coupling strength g. First of all, 
we choose the frequencies u\ and w 2 of the driving fields such that the detunings A 1 and A 2 
satisfy the resonance conditions: 

Ai =e liP -e 0>p = -g 2 /ui M (7) 
A2 = — £o,p+i = — g 2 /wm (8) 

In this way the driving field fli can resonantly couple the states \ipo,p) an d as illustrated 
by vertical arrows in Fig. lb. Similarly, the driving field Q2 can resonantly couple the states 
\ip 0>p+ i) and \ipx iP ). 

It is important to note that e nyP depends nonlinearly on the photon number n, therefore 
the above resonance relations do not hold for states with photon numbers larger than 1. In 
fact, if the optomechanical coupling strength g is sufficiently strong, the conditions (7) and 
(8) correspond to far off resonance for transitions from 1-photon states to 2-photon states, 



and this leads to the photon blockade effect 25]. Hence if the driving fields are sufficiently 



weak, it is justified to ignore the off-resonant transitions from 1-photon states to 2-photon 
states. Specifically, we require (for i = 1,2), 

ni4Z.\2g 2 /u M -Ku) M \ (9) 

where K is the nearest integer to 2(#/u;m) 2 - For example, if g < then K = 0. 

The inequality (9) means that f2j should be much less than the smallest detuning between 
2-photon manifold to 1-photon manifold. 

With the conditions (7-9), if the initial cavity photon number is zero, the system can be 
effectively confined to the subspace of zero and one photon. Furthermore, we will keep only 
the co-rotating terms in (4) as a rotating wave approximation. Under such an approximation, 
Qi can drive the corresponding resonant transitions only, and the Hamiltonian (4) in the 
interaction picture [i.e., the first term in (4) is removed by a rotating frame] becomes, 

N-l 

K = Y1 [{ A pI^ \M (M + IVO (^o, P+ i|) + h.c] . (10) 

Here the upper limit of phonon number iV can be infinite in general. However, a finite value 
of TV is possible if the transition between |V>o,tv) an d |^i,jv) can be completely suppressed 
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FIG. 2: (Color online) Solution of satisfying Eq. (11) as a function of N. For each N, there 
are multiple roots, and only the smallest positive root is shown. 

(Fig. lb). This can be achieved by a specific value of g = such that 

A ( i] N = eM-A)L% [qM) = 0. (11) 

In other words, one can design the Hamiltonian H' r with a prescribed N by using the 
optomechanical coupling strength g^. In Fig. 2, we illustrate the roots g^ as a function of 
N. We see that the value of g n decreases when N becomes larger and larger. For example, 
g N « 0.12w M for = 100. 




IV. DARK STATES 

The Hamiltonian H' r has an eigenvector \D) with a zero eigenvalue: 



N 

\D) = CY,Mp)m®\0) c (12) 

p=0 



where {3 = 1, and for p > 0, 



n 2 J 11 a (1) 

and C is a normalization constant. Such an eigenvector is a coherent superposition of phonon 
states and the cavity contains no photon. It is a dark state induced by quantum coherence 
of the mirror, and the interference forbids any excitation of cavity field even though the 
cavity is constantly driven by the two external fields. 
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FIG. 3: (Color online) Phonon number distributions of dark states for different ratios of driving 
strengths at the coupling strength g = g\Q = 0.37ujm- 
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FIG. 4: (Color online) {(An) 2 )/(n) as a function of fii/0 2 for g/u M = 0.17, 0.64, 0.76 corre- 
sponding to N = 20, 3, 2 respectively. 

The phonon number distribution of the dark states is complicated by the Laguerre func- 
tions in Eq. (13). While the details form of \(3 p \ 2 requires a numerical evaluation of Eq. 
(13), we find that \(3 P \ 2 is mainly controlled by the ratio of the strengths of driving fields. 
Such a feature is illustrated in Fig. 3 for various f^/^i- For example, when f^/^i = 3, the 
probability decreases quickly with the increase of phonon number m. In the case f^/ili = 1, 
there is a peak appears in the probability distribution. If the ratio is further decreased to 
^2/^1 — 1/3 , the peak is shifted towards higher phonon numbers. 

It is worth noting that phonon number distributions of dark states exhibit a sub- 
Poissonian statistics. This is illustrated in Fig. 4 in which the ratio ((An) 2 )/(n) as a 
function of Q1/Q2 is plotted. We see that ((An) 2 ) / {n) is always less than 1 (i.e., sub-Poisson 
distribution) except at the small region Qi/f^ near zero. In particular, the ((An) 2 )/(n) de- 
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creases with Oi / 2 . In Fig. 4, we also find that the curves are quite insensitive to the value 
of gN used. 

V. PREPARATION OF DARK STATES BY CAVITY-FIELD DAMPING 

In this section we discuss how the system can be optically pumped into the dark state 
by cavity-field damping. The evolution of the dissipative system is governed by the master 
equation: 

^j- = -i [H r , p] - — (a f ap - 2apa) + pa ] a) (14) 

where p is the density matrix of the photon-mirror system, and 7 C is the cavity-field decay 
rate. Here we have assumed that the mechanical motion of the mirror has the damping 
rate j M that is much smaller than j c . In addition, we will focus on a finite time interval 
l/7c «C t <C 1/7m during which the optically pumping is complete while the decay of mirror's 
motion remains negligible, and so we only include the decay effect of the cavity field in the 
above master equation. Note that the original Hamiltonian H r defined in Eq. (2) is used 
without employing the rotating wave approximation in Eq. (10). If H r is simply replaced 
by H' r , then p = \D)(D\ is already a steady state solution of the master equation Eq. (14), 
because the cavity field damping term has no effect on \D) (which has zero photon). 

We have solved the master equation (14) numerically with an initial ground state of the 
system. Specifically, we are interested in the fidelity F defined by 

F = Tr(\D)(D\p(t)), (15) 

which measures the probability of the system in the dark state \D). Some examples are 
given in Fig. 5 in which F increases with time and approaches a steady value close to 1 
in a finite time. For the three cases shown in Fig. 5, the fidelities can reach F 0.99 
with g N /u M = 0.37 (JV = 10) at the time T w 8000^/. To ensure that the mechanical 
decoherence is negligible, we need 7m < 10~ 4 om for the parameters used in Fig. 5. Indeed, 
we have tested numerically the performance by including mechanical damping in the master 
equation, and found that F 0.99 when 7m = 10~ 5 um and F ps 0.93 for a larger 7m = 
10~ 4 Co>m- We remark that the time T required to generate dark states is shorter for smaller 
N. For example, when N = 3, our numerical calculations indicate that T ps 2000u; m 1 . 
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FIG. 5: (Color online) Time evolution of the fidelity F for various O2/O1 ratios. The parameters 
are: N = 10, g/oj M = 0.37, j c /co M = 0.05, Q 2 /u M = 0.01, Ai/wju = -0.14, A 2 /w M = -1.14. 

The increase of F is understood as an optical pumping effect. This is because when a 
photon leaks out of the cavity, the mirror can have a non-zero probability making a transition 
to the dark state. Since the dark state is decoupled from the driving fields, it can no longer 
be excited, and hence the occupation of the dark state accumulates as time increases. In 
our system, the loss is mainly due to a leakage of phonon population beyond the phonon 
number N, because the cavity field decay also causes the mirror to make transitions to states 
other than the dark state. Note that the transition amplitudes from |1) c |p)m to \0) c \q)M 
due to the transmission of a photon out of the cavity is proportional to the Frank-Condon 
factor (p\q)M, the loss can be reduced by choosing a sufficiently high N and Q 2 > Qi. This 
is because dark states with Q 2 > ^1 concentrate on lower phonon numbers (Fig. 3) and 
hence the Frank-Condon factors for transitions to states of phonon number higher than N 
is smaller. 

VI. CONCLUSION REMARKS 

To conclude, we have addressed a quantum interference effect in the single-photon strong- 
coupling regime of optomechanics. In such a regime we discover a class of dark states of a 
moving mirror under the conditions (7-9). These dark states make the cavity decoupled from 
two external driving fields, and the decoupling is derived without employing the linearization 
treatment. We provide an analytical expression of the dark states which indicate the depen- 
dence of the ratio Q 2 /Q\ of the driving fields and the optomechanical coupling strength g. 
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With the help of cavity damping, dark states can be prepared by optical pumping. 

In this paper specific optomechanical coupling strengths have been used in order to 
'trap' the mirror state in the space of finite phonon numbers. It is important to ask how 
a slight deviation of g away from g^ would affect the dark state generation, because the 
trapping effect would become imperfect when g ^ g^. To address the issue, we have tested 
the sensitivity of the fidelity to small variations of g values. For example, with ^2 > ^1 and 
parameters in Fig. 5, by introducing a 3% deviation of the g value we found that the fidelities 
can still reach about F ~ 0.99. This can be understood by Eq. (13). Strictly speaking, 
when g is not exactly equal to gN, the upper limit N in Eq. (12) should be extended to 
infinity because A$ N is no longer zero. However, if f^/f^ is less than one, then (5 p appearing 
in Eq. (13) can quickly converge to zero, yielding a significant population concentrated at 
phonon numbers much lower than iV (Fig. 3). In this case the leakage of phonon numbers 
beyond N would be negligible and hence a slight deviation of g would not affect the dark 
state generation appreciably. On the other hand, we find that the system with ^2 < ^i is 
quite sensitive to g. This is because if Oi/f^ > 1, P P in Eq. (13) tends to shift to higher 
phonon numbers (Fig. 3), making the leakage beyond N significant. Therefore it would be 
more favorable to employ Q2 > ^1 for experimental realizations of our scheme. Finally, we 
remark that the realization of single-photon strong-coupling regime is still a challenge for 
current experiments, but some recent progress have been made in various optomechanical 



systems in order to reach the strong coupling regime 29j-|3j|. 
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